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We study secondary electron emission from metallic surfaces due to Auger de-excitation of di- 
atomic metastable molecules. Our approach is based on an effective model for the two active 
electrons involved in the process - a molecular electron described by a linear combination of atomic 
orbitals when it is bound and a two-center Coulomb wave when it is not and a metal electron de- 
scribed by the eigenfunctions of a step potential - and employs Keldysh Green's functions. Solving 
the Dyson equation for the retarded Green's function by exponential resummation we are able to 
treat time-nonlocal self-energies and to avoid the wide-band approximation. Results are presented 
for the de-excitation of N2( S«) on aluminum and tungsten and discussed in view of previous exper- 
imental and theoretical investigations. We find quantitative agreement with experimental data for 
tungsten indicating that the effective model captures the physics of the process quite well. For alu- 
minum we predict secondary electron emission due to Auger de-excitation to be one to two orders of 
magnitude smaller than the one found for resonant charge-transfer and subsequent auto-detachment. 

PACS numbers: 34.35.+a, 79.20.Rf, 79.20.Hx 



I. INTRODUCTION 

De-excitation of metastable atoms and molecules with 
simultaneous release of an electron is a surface scat- 
tering process of great technological importance. De- 
excitation of atoms is used as a surface-sensitive elec- 
tron spectroscopy-^ - — and de-excitation of molecules is 
an important process in molecular low-temperature gas 
discharges. It is one of the main wall-based secondary 
electron emission channels controlling, together with wall 
recombination and various volume-based charge produc- 
tion and destruction channels, the overall charge balance 
in the discharged In the de-excitation process both the 
target and the projectile are composite objects. A great 
variety of reaction channels is thus conceivable making 
the investigation of this scattering process a challenging 
t ask, particularl y for molecules. 

IStracke et alJ & experimentally investigated the de- 
excitation of metastable nitrogen A^( 3 £^) molecules on 
a tungsten surface and proposed two main reaction chan- 
nels. Firstly, the Auger de-excitation (also referred to as 
Penning de-excitation) , 

7V 2 ( 3 E^ + e „ l ^iV 2 ( 1 S+)+ e/ , (1) 

where e m and e / denote an electron inside the metal and 
a free electron, respectively, and secondly, the formation 
of the A^( 2 n 9 ) shape resonance with subsequent auto- 
detachment, 

N 2 ( 3 Xt) + e m -> A^( 2 n 9 ) -> iV 2 ( 1 S+) + e f . (2) 

IStracke et alj £ concluded that out of these two compet- 
ing processes reaction ([2]) should be more efficient, as 
it is a combination of two single-electron charge-transfer 
transitions, whereas ([T]) represents a less probable two- 
electron transition. Using thermal molecules they mea- 
sured the energy spectrum of the released electron and 
estimated the overall secondary electron emission coeffi- 
cient 7 e , that is, the probability for releasing an electron 



by a single metastable molecule de-excitation at the sur- 
face, to be about 1CU 3 — 1CU 2 . The experimental esti- 
mate for 7 e does not discriminate between the two reac- 
tion channels . It rather includes both channels. Indeed, 
IStracke et alj ^ mention that in the spectrum of the emit- 
ted electron they also observe a weak signal due to Auger 
de-excitation. It is one order of magnitude weaker than 
the signal due to charge-transfer. 

Based on the assumption t hat the charg e-transfer 
channel ([2]) is the dominant one lLorente et alj 7 - theoreti- 
cally investigated the de-excitation of A^ 3 !^) molecules 
on an aluminum surface. The resonance-driven sec- 
ondary electron emission coefficient resulting from their 
calculated electron emission spectrum is about 1CU 1 
which is one order of magnitude larger than the value 
IStracke et alj £ give for tungsten. The Penning chan- 
nel (UJ was not considered by Lorente and coworkers. 7 
Its strength for an aluminum surface is thus unknown. 

In the pres ent work we a dopt the point of view com- 
plementary to lLorente et al.l 's— investigating Penning de- 
excitation while neglecting any contribution due to res- 
onant charge-transfer. In particular for tungsten the ef- 
ficiency of the Penning process may be comparable to 
the efficiency of the charge-transfer process, because the 
molecular orbital hosting the hole in the electronic config- 
uration of Af 2 ( 3 E+) is roughly 2.5 eV below the bottom of 
the conduction band of tungsten^ To bring this orbital 
in resonance with conduction band states of the metal re- 
quires therefore a large image shift and broadening due to 
the interaction with the metal. Rough estimates of these 
two effects based on what is known about them for alkali 
atoms interacting with surfaces^ imply that the reso- 
nance condition can only be met for vibrationally excited 
A r 2 ( 3 S+) states. Thus, at least for A^ 2 ( 3 S+) in its vibra- 
tional ground state, Penning de-excitation and charge- 
transfer are eye-to-eye competitors. For aluminum the 
electronic band structure is much more favorable for the 
charge-transfer scenario.™ Here, the bottom of the con- 
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duction band is only 1 eV above the molecular orbital in 
question. This energy difference may be bridged by the 
combined action of image shift and level broadening. 

In order to theoretically analyze the Auger de- 
excitation of diatomic molecules on metallic surfaces we 
adopt an effective two-electron model, where one electron 
resides in the conduction band of the metal and the other 
in the excited state of the molecule. The metal is mod- 
eled as a half-space containing a free electron gas charac- 
terized by a work function and a conduction band depth, 
while the molecule is modeled in terms of a two-level sys- 
tem with energy spacing corresponding to the excitation 
energy of the metastable state. Focusing on the situa- 
tion where the emitted electron (Auger electron) origi- 
nates from the molecule we describe it by a two-center 
Coulomb wave. The coupling between the molecule and 
the metal is through the Coulomb interaction between 
a metal electron and the electron in the excited state of 
the molecule. Image interactions are neglected (except 
for the emitted electron, where we include it in terms of 
a surface transmission function) because they give rise to 
a hybridization of metallic and molecular single-electron 
states^ which we prefer to discuss in connection with 
direct charge-transfer. 

Although our model for Auger de-excitation is some- 
what crude, it contains the most relevant degrees of free- 
dom in a reasonable approximation and can be param- 
eterized by energies which are relatively easy to obtain. 
In particular the latter aspect is rather important for us, 
because our interest in Auger de-excitation and related 
processes stems from their relevance for bounded low- 
temperature gas discharges. It is well known that elec- 
trons can be released from plasma walls by de-excitation 
of metastable species but for most wall materials and 
projectiles the secondary electron emission coefficient is 
unknown. A flexible, easy-to-use microscopic model for 
its calculation is thus needed. 

Following the lead of Makoshi and coworkers ; 12 ' 13 who 
investigated the de-excitation of metastable atoms, we 
employ in the following the Keldysh technique^— to 
calculate within the trajectory approximation^ the sec- 
ondary electron emission coefficient and the spectrum of 
the emitted electron for a diatomic metastable molecule 
hitting a metallic surface. A description of this type 
of surface collision with Green's function s 12 ' 13 ' 19 " — is 
mathematically more demanding than using rate equa- 
tions j 18 ' 23 ~ — Green's functions are however rather flex- 
ible in handling the non-adiabaticity of the projectile's 
motion; 12 ' 13 the Coulomb correlations on the projec- 
tile ji2r— and the collective electronic excitations of the 
surfaced In addition, image shifts and level broaden- 
ing due to image interactions as well as vibrations of the 
molecule may also be included in a theoretical description 
based on Green's functions. 

In contrast to Makoshi's work ) 12 ' 13 our approach is not 
restricted to time-local Auger self-energies and thus to 
the wide-band approximation. To overcome this limita- 
tion we solve the Dyson equation for the retarded Green's 



function by exponential resummation^ Our approach is 
also not restricted to phenomenological Auger matrix el- 
ements. We work with the full matrix element, exploit- 
ing only the locality of the bound molecular states and 
the large distance of the molecule's turning point from 
the surface. Although the final equations for the sec- 
ondary electron emission coefficient and the spectrum of 
the Auger electron are highly complex they can be nu- 
merically evaluated within an interpolatory grid-based 
Monte-Carlo integration scheme. 

We specifically apply our approach to Auger de- 
excitation of A2( 3 £^) on an aluminum or a tungsten 
surface. The metastable molecule is assumed to be in 
its vibrational ground state and the molecule's turning 
point is obtained from the surface potential applicable to 
A2( 3 S+) on metallic surfaces. For an aluminum surface 
we find the secondary electron emission coefficient due to 
Auger de-excitation ([T]) to be one to two orders of magn i- 
tude smaller than the one deduced from lLorente et alj 's^ 
theoretical study of the direct charge-transfer channel © 
and for tungsten we find good quantit ative agreeme nt 
with an experimental estimate based on lStracke et alj 's^ 
measurements. 

The paper is structured in the following manner. In 
the next section we give more details about the effec- 
tive model on which our analysis of Auger de-excitation 
is based. In Sect. [HI] we employ the Keldysh technique 
to extract physical quantities from the model. Thereafter 
we describe in Sec. II VI the numerical scheme we employed 
for the calculation of the secondary electron emission co- 
efficient and the energy spectrum of the Auger electron. 
Results are presented in Sec. |V] We conclude the pa- 
per in Sec. IVII and complement it by two appendices. 
Appendix contains the explicit form of the wave func- 
tions we used in our calculations and Appendix [B] fixes 
the notations of the Keldysh formalism. 



II. MODEL 

We investigate the de-excitation of a metastable nitro- 
gen molecule impacting on a metallic surface with simul- 
taneous release of an electron. The model we employ is 
an effective one that concentrates on the most important 
degrees of freedom and enables us to describe the system 
by a few parameters which arc accessible through exper- 
iments or theoretical calculations. The primary goal will 
be to calculate the secondary electron emission coeffi- 
cient 7 e . 

Focusing on the essentials of the process, we intro- 
duce from the start some simplifications and restrictions. 
First, we assume the metal surface to be planar, ideal, 
uncharged, and to stretch over the entire half space z < 0. 
Furthermore, we consider only the dominant metastable 
state -Zv2( 3 Ej). In addition, we employ the trajectory 
approximation^ that is, we decouple the translational 
motion of the molecule from the dynamics of the sys- 
tem and externally supply its trajectory. Finally, the 
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FIG. 1: Schematic illustration of the collision geometry. 

molecule is assumed to impact the surface under normal 
incidence with constant velocity v and constant angle tp 
of its axis to the surface. Because of the translational 
symmetry of the solid surface in the x-y-plane, it is then 
sufficient to consider only rotations of the molecule axis 
about one particular axis in this plane, for instance, the 
y-axis (see Fig[T|). 

We now cast the model assumptions into mathemat- 
ical form, starting with the trajectory. Assuming the 
molecule to start moving at to = —oo and to hit the 
surface at t = the trajectory of its center of mass is 

R(t) = z R (t)e z = (v\t\+z )e z , (3) 

where zq denotes the turning point. The center of mass 
motion is classical. Hence, the turning point zq can be 
determined by considering the motion of the molecule in 
the molecule-surface interaction potential Vs(z) for given 
initial kinetic energy Ekin- Using a Morse- type potential, 
energy conservation gives for the position of the turning 
point 

*'=-~<- l A i+ f^f) w 

with material parameters d, a, and z e .— For metals 

d = 0.738 eV, a = 640pm" 1 , z e = 245pm. (5) 

According to the authors of Ref. [2^ these values are not 
very specific to the particular metal. 

To set-up, for a given trajectory, a Hamiltonian for 
a nitr ogen molecule de-exciting at a metal surface we 
follow ICazalilla et alJ 2£ and distinguish between indirect 
de-excitation (Penning process, solid lines in Fig. [2]), in 
which the electron is emitted from the molecule, and 
direct de-excitation (exchange process, dashed lines in 
Fig. [2]), in which the electron is emitted from the surface. 
It must be stressed that although the asymptotic form 
of the wave function of the emitted electron is in both 
cases a plane wave, in the region relevant for the calcula- 
tion of the matrix element, the Auger wave function for 
the Penning process resembles a single-electron contin- 
uum wave function of the nitrogen molecule whereas the 
Auger wave function for the exchange process is basically 



a continuum state of the solid with positive energy. In 
the following we will only consider the Penning process 
because our calculations showed that its matrix element 
is much larger than the matrix element associated with 
the exchange process, a consequence of the orthogonality 
of the bound molecular wave functions. 

We construct the Hamiltonian by combining three dif- 
ferent kinds of single-electron states to a single-electron 
basis: the single-electron states of the conduction band 
of the solid surface, which we approximate by the states 
corresponding to an electron trapped by a step potential 
of depth $cr the free single-electron states associated 
with the molecule's continuum for which we use a two- 
center Coulomb wave^^— and effective single-electron 
states for the bound states of the molecule. To keep the 
description of the molecule as simple as possible we ap- 
proximate the latter by a degenerate two-level system 
keeping, within the LCAO representation of the nitro- 
gen molecule* 3 - 3 - only the 2ir u and the 2ir g molecular or- 
bitals (MOs) which are the two MOs whose occupan- 
cies change during the de-excitation process. To con- 
struct the two MOs we use moreover hydrogen-like wave 
functions with effective charges to mimic the Roothaan- 
Hartree-Fock wave functions of atomic nitrogen* 3 - 4 In the 
molecule's ground state A"2( 1 E+) the 2ir u MO is fully oc- 
cupied and the 2ir g MO is empty while in the excited 
state 7V2( 3 S+) the 2ir u MO contains a hole and the 2ir g 
MO is singly occupied. Both of these levels can carry four 
electrons and are degenerate in the electron spin s = ±i 
and the magnetic quantum number m = ±1. Since the 
process we consider does not involve any spin flip, we 
ignore the spin. We can thus label the ground state of 
the two-level system and its excited state by 0m and lm, 
respectively, and denote the corresponding energies by 
eo and e\. The states of the metal and the free states 
are labelled by k and q, respectively. The mathematical 
expressions for the wave functions of the single-electron 
states are given in Appendix A. 

The description of the electronic structure of the 
molecule-surface system is completed by aligning the 
single-electron states against each other and against the 
vacuum level by use of the metal's work function 
the metal's conduction band depth <&c, the molecule's 
ionization energy Ae^, and the excitation energy of the 
molecule Ae e . The metal states are of course occupied up 
to the Fermi level ep. Our model is thus characterized 
by a few energy parameters, an effective charge, and a 
bond length (which enters the molecular wave functions 
of the diatomic molecule) . 

The electronic structure of the simplified model is 
sketched in Fig. [2l together with the transitions of the 
Penning de-excitation and its associated exchange pro- 
cess. Due to the symmetries of the molecular ground 
state ( 1 S+) and the molecular excited state ( 3 S^), only 
transitions with Am = are involved. They are driven 
by the Coulomb interaction between the excited electron 
in the 2n g MO and an electron in the Fermi sea of the 
metal (see, for instance, Refs. [24Tj27l ). The three elec- 
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FIG. 2: Energy scheme of the simplified model showing 
Penning de-excitation (solid lines) and its exchange process 
(dashed lines). Also indicated are the classical energy cut- 
offs £g" n and e^ ax which can be calculated from the energy 
balance ei + eg = eo + £q that holds in the adiabatic limit. 



trons in the 2tt u MO act only as spectators and can thus 
be neglected. Assuming moreover the Fermi surface of 
the metal to be rigid, the de-excitation of A^( 3 £j) is 
basically a two-body scattering process, whose Hamilto- 
nian, written in the single-electron basis described in the 
previous paragraph, is given by 



H=H + Hi(t), 

k 9 

rn m 

Hi(t) = £ 0(*)4mC fe -4c lm + H, 

k,q.7n 



(6a) 
(6b) 

(6c) 



where Hq represents the Hamiltonian of the non- 
interacting system, and Hi(t) contains the Penning de- 
excitation. 

The Auger matrix element V n 9,1 '? contains the time- 

° m,, re 

dependence of the Hamiltonian, and thus carries the in- 
trinsic non-equilibrium character of the system. In terms 
of the single-electron states given in Appendix A, it can 
be written as 

W) = J&Jd? *U*V(*))*#) (7) 
x V c (\r- r'|)*t (^(t))¥ lm (^(t)) , 

where Vc represents the Coulomb potential and the sub- 
script ip denotes the associated vector as seen from the 
molecule's reference frame, which is centered about the 



molecule's center of mass R(t) and has its z-axis aligned 

along the molecule axis. The vectors f$ (i) and q v are 
thus given by 



q v = Cl((p)q , 



(8a) 
(8b) 



where the matrix Q((p) describes the rotation around the 
y-axis (see Fig. [T]). Due to the diatomicity of the nitrogen 
molecule the interaction matrix element depends on the 
orientation of the molecule with respect to the surface, 
that is, on the angle ip. For convenience we suppress 
however this dependence in our notation of the matrix 
element. Note, since we are considering the Penning pro- 
cess, the continuum states q are defined with respect to 
the molecule's center of mass. Thus the associated wave 
function within Eq. (J7J must depend on the molecular 
reference frame variables q v and r'(t). 

As already mentioned, we take the image interaction 
only for the emitted electron into account. It always 
feels the image potential and can thus only escape from 
the surface when its perpendicular energy is sufficiently 
highii Therefore, in the calculation of the secondary elec- 
tron emission coefficient and the spectrum of the Auger 
electron we will multiply the matrix element by a step 
function 



e 
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167T£ Z R (t) - Z t 



(9) 



where Zi is the position of the image plane. For alu- 
minum and tungsten Z{ is given in Table U and turns 
out to be closer to the surface than the turning point of 
the molecule. Equation ([9]) is thus well-defined along the 
trajectory of the molecule. It is sometimes referred to as 
surface transmission function,— although for the case of 
Penning de-excitation this term is somewhat misleading, 
as the escaping electron is not emitted from the inside of 
the solid but from the molecule. 



III. QUANTUM KINETICS 

The model established in the previous section will now 
be treated using the Keldysh technique, a brief descrip- 
tion of which is given in Appendix FBI 

We start by calculating the unperturbed Green's func- 
tions Q^a , with a and /3 representing any of the labels we 
used to characterize the single-electron states of the sys- 
tem, k, q, Om, or lm. Since the time evolution of the free 
Green's functions is determined by Hq, they are diago- 
nal, that is, Q^l ~ 5 a p. For convenience we abbreviate 
the double subscript aa by a in the following. Inserting 
the solutions of the interaction-free Heisenberg equations 
for the creation and annihilation operators appearing in 



5 



model ([6]) , the free propagators read 

iG* w {t,t') = e(i-t')e* EQ( ' _t ' ) , (10a) 
iG^\t,t') = -e(t' -t)^ 5 ^-^, (10b) 
iG^°\t,t') = [1 - 2n a (to)] e* 6 ^*-*'), (10c) 

where n a (to) is the initial occupancy of the state a at 
to = — oo. 

In accordance with the model we introduced in the pre- 
vious paragraph we assume the excited molecular level 
to be initially occupied with a single electron of mag- 
netic quantum number m — fi, that is, ni m (to) = S rnfl . 
The molecular ground state level, in contrast, is empty 
at to because we neglect the spectator electrons. Hence, 
nom(to) = for all m. The free electron states are also 
empty at to, implying n${to) = 0, and the electronic 
states within the metal are initially filled up to the Fermi 
energy sp, that is, ng(to) = @(ef — eg). 

For the calculation of the full Green's functions G a /3 we 
need expressions for the se lf-energies £ aj g, which, in line 
with the work by iMakoshJ ^ we derive from a diagram- 
matic expansion up to second order in the Auger matrix 
element. This is justified because the Auger matrix ele- 
ment is in general very small. Because of the diagonal- 
ity of the unperturbed Green's functions the self-energies 
and thus the full Green's function are also diagonal. 

We first investigate the excited molecular state. Fig- 
ure [3] shows the only non-vanishing second-order self- 
energy diagram for Ei M . It can be evaluated to 



Ef /t (t 1 ,t 9 ) = e(ti-t9)sf TO (ti,*9), 

S^(*i,ta)= -e(* 2 -ti)EL(*i>«2) 



(lla) 
(lib) 



= - «e raw]* v ::m 

(He) 

x n j :(to)e-i ( - eo+ ^- £k)( - tl - t2) . 



Using Eq. (|lla|) . the Dyson equation for the retarded 
Green's function (Eq. (IBlla[) in Appendix B) can be 
solved iteratively. The result is 



G^(t,t') = Gf°\t,t')W^t,t') 



with the infinite series 



w fl (t ) t , ) = y^w^(t,f) 



!/=0 

whose individual terms, , are given by 
ft rt 



(12) 



(13) 



W^{t,t') = {-!)" / dtj dt 



2 • ■ • / dt 2v 
v Jv 



v jv jv (14) 

x A |t (ti,t 2 ) . . . A^(t 2u -i,t 2 u) , 

where we introduced the quantity 

A„(t 1 ,t 2 )=iXK(t 1 ,t 2 )e^-^ , (15) 




FIG. 3: Diagrammatic representation of the self-energy 
tz) of the excited molecular state in second order per- 
turbation theory. 



which emerges from the self-energy terms of the iterated 
Dyson equation. 

The infinite series (fT5)) is exact but useless. To obtain 
an expression for the retarded Green's function which 
is amenable to further manipulations we employ the 
exponential resummation technique (see, for instance, 
Ref. [H). For that purpose we introduce a new func- 
tion F/j,(t,t') and perform a perturbation expansion of 
and in terms of A M . Using the virtual expansion 
parameter A = 1 we write 



e ^n = e ^ F ^\ (16) 



!/=0 



Expanding the exponential in (|16[) and then comparing 
the different orders of A leads to explicit expressions for 
the expansion coefficients Fa , the first few of which are 



= W (1) 

F {2) = w {2) 



(17a) 
(17b) 

(17c) 



where, for convenience, we dropped the time arguments. 

The retarded Green's function can now be conveniently 
written as 



G«(t,t>) = G%°Xt,t>)eW') 



(18) 



with the function F^ in the exponent given by the sum 
of the terms in (fl"7|) . Using relation (|B10[) together with 
Eq. (fl"8|) the advanced Green's function becomes 



Gt(t,t') = G^ \t,t')ei^'^ 



(19) 



To calculate the occupation of the excited molecular 
state we also need the Keldysh part of the Green's func- 
tion defined in Eq. (|B12[) of Appendix B. Using the ex- 
plicit form of the free Green's functions (fit))) we first 
rewrite this equation into^ 



G£(f,t / ) = -i [1 - 2n a (to)]G%(t,to)G£(to,t') 
+ [ ' dt x fdt 2 G B a (tM)^{hM)Gt{t 2 ,t') , 

J t Jt 



(20) 
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using the identity 
G* (0) (M') 



-»[l-2n„(to)] 
xGf)(Mo)Gf)(i„,i'), 



(21) 



and the Dyson equation of the retarded and advanced 
Green's function (|Bllal) . Note that Eq. (J20J) is not 
limited to the excited molecular level. It holds for all 
states a. 

Inserting (jllcl) . (fig]) , and dTHJ) into (gO]) we can now cal- 
culate the Keldysh part of the Green's function G^(t, t'). 
Taking the latter at equal times t = t' and utilizing 
Eq. (|B13|) we finally obtain for the occupancy of the ex- 
cited level at time t, 



(t) 



= 2SR[F M (t,* )] 



(22) 



where 3?[. . . ] denotes the real part. To lowest order in the 
interaction, that is, to lowest order in the occupation 
of the molecular excited state is given by 

n<J(i) = e 2JW <*.*>)] = e-J^f&tz _ (23) 

We now turn to the free electron states, that is, the 
states which may get occupied by the electron released 
by the de-excitation of the molecule. A treatment of 
these states analogous to the excited state leads to the 
following expression for the occupation at time t 



(24) 



where is defined the same way as F^ but with A M 
replaced by A M ^. The latter is implicitly defined through 
the relation 



(25) 



using Eqs. (J3]) and ([Tic]) for A fl (t 1 ,t 2 ). 

Equation (j2~4")l is not very useful, because it can- 
not easily be summed over q, which is however needed 
to calculate the secondary electron emission coefficient 
7 e . Because of this obstacle we adopt the approach 
of Makosh i 12 ' 13 and go back to Eq. (f20|) . expand it for 
equal times t = t' and a = q up to first order in the 
self-energies, and insert the result into Eq. (|B13|) which 
yields^ 3 - 



ntf) = -i [ dh f dt 2 Gf°\t,t x )^{txM) 

Jto Jto 

xG?°\t 2 ,t>). 



(26) 



To account in Eq. (|26|) for lifetime effects o f the 
metastable molecule we follow once more iMakoshj ^ and 
employ in the calculation of the self-energy the full 

(" dressed" ) Green's function of the excited state Gf~ in- 
stead of the unperturbed one. Up to second order (see 




FIG. 4: Diagrammatic representation of the self- 
energy T,^(t\,t2) of the emitted electron in second order 
perturbation theory. The dressed Green's function Gi M is 
indicated by a double line. 



Fig. 0} we obtain 



x G 0A1 +(0) (t 2 , tl ) Gt- (0) (H , t 2 ) G+- ( tl ,t 2 ). 



(27) 



In order to proceed we need to calculate G^T . For that 
purpose, we first transform the corresponding component 
of the matrix Dyson equation (IB7[) into 



G 



G 



-(0) 



Gfr^Gt-, (28) 



employing the fact that up to second perturbation order 
= Ef M and T,^~ = — . Inserting the free Green's 
functions (|10[) as well as the self-energies (fTT]) and the full 
advanced Green's function (|19p one can solve Eq. (j2"5)l 
itcratively to obtain 

G+-{t,t') = G+- [0 \t,t')e F ^^ e\ F ^'^Y . (29) 



After inserting Eq. (Qjj and ([29]) into Eq. (|27j) we obtain 
for the Keldysh part of the self-energy 



E f (*i.*a)=^E V o^(*0[ V o^(*a)]" 



x nAt ) e *( e o-e*-ei)(ti-t a ) 



with 



V f,1 i*ft) = V q '^(t) 



(30) 



(31) 



the renormalized Auger matrix element.— 

To compute finally the occupation of the g-states, we 
insert Eq. ([30f into Eq. (|26|) and obtain 



(32) 



dh f dt 2 A^(tx,t 2 ) 



to 



where A M ^ is defined in Eq. (|25[) with the plain matrix 
elements replaced by the renormalized matrix elements 
leading to 



A„g-(ti,t 2 ) = A„g-(ti,i 2 )e 



[^(ti.to)]*^^,^) 



(33) 
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Equation (1321) represents the spectrum of the emitted 
electrons. The secondary electron emission coefficient 7 e 
can be calculated from Eq. (|32l) by taking t = oo and 
summing over all possible q, 

/•oo />oo 

7e = V%(oo) = / dt x / dt 2 A p (ti,i a ) (34) 

- ./to Jto 

with A M defined by 

A M (ti,t 2 ) = 53^.*(* 1 '* 2 ) ■ ( 35 ) 

Equation (|34[) can be rewritten to 

Te = 1 - e 2^(oo,to)] . (36) 

The lowest two orders of Eq. f[36j) in terms read 

/•OO /> CO 

7 (°) = dti dt 2 A^(ti,t 2 ) , (37a) 
7 (i) =1 _ e-St&Stp* ^^A) . (37b) 

The preceding calculation of the self-energies does not 
treat free and excited states on an equal footing. Only the 
self-energy for the free states is renormalized whereas the 
one for the excited state is not. As a result, particle con- 
servation is not strictly guaranteed when the correspond- 
ing occupation n umbers a re calculated. The same short- 
coming holds for lMakoshil 's^ original approach. Our nu- 
merical results showed however that particle conservation 
is obeyed for all physically relevant situations, justifying 
the treatment of the self-energies a posteriori. 

Let us finally remark that although, as far as the 
logic of our approach is concerned, we closely followed 
iMakoshil ^ our results have wider applicability. In con- 
trast to him we do not work with a real phenomenological 
Auger interaction, depending only on time, and do not 
employ the wide-band approximation for the free states 
which would lead to time-local self-energies. We are also 
not restricted to the lowest order expressions given in 
Eq. ([23]) and {(37]) . In principle, we can calculate correc- 
tions to these expressions using the higher order expan- 
sion coefficients of given by Eq. (TT7)) . 



IV. NUMERICAL SCHEME 

Inspecting the main results of the previous section, 
Eq. (f22|) and (f36|) . we realize that these equations are 
highly complex. For instance, calculating A p from 
Eq. (|I5p requires summation of a matrix element prod- 
uct over all k and q vectors, which, letting the box size 
L — > oo, equals a sixfold integral. The Auger matrix el- 
ement itself, according to Eq. ([7]). involves another six 
dimensional integration over f and r'. Since in Eq. (|I5p 
the matrix element occurs as a product at two different 
times, this makes a total of eighteen dimensions of inte- 
gration. Calculating the level occupancies from Eq. (|22|) 
and (|36l) requires at least another two-dimensional inte- 
gration over the time arguments of A^. Thus, employ- 
ing these equations as they stand for numerical calcula- 
tions is clearly out of reach if reasonable computing time 
and considerably small numerical errors are required. In 
this section we will therefore introduce an approximation, 
which makes those numerical calculations feasible. 

To simplify the matrix element (|7|) we first utilize the 
particular form of the wave functions (see Appendix [A]) . 
Since the molecular wave functions $ 0m / lm are localized 
on the molecule and the metal wave functions \I> > and 
the free wave functions \& g are bounded everywhere (in 
the mathematical sense), the main contribution to the 
matrix element ([7]) will arise from points close to the 
actual molecule position. Noting further, that for the 
low kinetic energies we are interested in (£/&„ < I eV) 
the molecule's turning point lies far outside the surface 
[zq > 4.35 as, see Eq. §4§ and the parameters ((5J). We 
can thus safely restrict the r-integration in (O to z > 0, 
that is, neglect the overlap of the wave functions inside 
the solid. Using then the transformation 

r^ = rt(<p)r m + R(t) , (38) 

along with the explicit form of the electronic wave func- 
tions (see Appendix [AJ we arrive at 

V£jt(t) « Cvn.e-^VVS , (39) 

with 



\N%Y 

C V = 1 8 — . ( 40a ) 



c 2 



47T£n 



V S =Jdn Jdf 2 Q{ziv + z ) ^^ ( e -l* + l +e -l*-l) (e-l^l - e -l*-l) ^ 
x e i(fc,*i„+*,vi*) e -«»„»i„ e -i&-ft C(q V) f 2+ ,Z)C(q (p ,f 2 -,Z) , 



8 



Material 


<E>w- [eV] 


$c [eV] 


Ref. 


Zi [as] 


Ref. 


Al 


4.25 


16.5 


2 


2.95 


35 


W 


4.5 


10.9 


jI 


3.0 


35 



TABLE I: Model parameters for aluminum and tungsten. 

where Qi and ipi (i = 1,2) are the usual cylindrical 
coordinates and we have introduced the abbreviations 
Z = Z c /2, r l± = n ± (5/2)e z (i = 1,2), and r ltp = 

(xi<p, Uiip, z \ip) — "HyOri- Here 5 is the bond length 
of the molecule. Note that in Eq. (|40bp we have used 
as the unit of length. Moreover, we did approximate 
Q{zi v + zn(t)) by Q{z\ v + zq), which, as numerical tests 
confirmed, is a good approximation for the distant turn- 



ing point positions we are considering. 

To calculate the occupation numbers from the equa- 
tions given in Sec. lIIIl we need to compute the time inte- 
gral of A fJi (ti,t 2 ). Since in Eq. (l39l) the time dependence 
has been separated we can insert the trajectory Q and 
carry out the time integration analytically. For the lowest 
order equations it suffices to calculate the integral with 
equal upper boundaries, that is, 

T(t) = fdh fdt 2 A p (ti,t a ) (41) 

Jto Jto 

Using spherical coordinates for the wave vectors k and q 
we obtain 



2(t) = 



9 (0 ^ M M H 2 V 2 k 2 £ 2 I ^ ^ k / ^ / ^ / M 

2(27r) A , 2 7r „A , 2 7r " v K v £ o Jo Jo Jo Jo Jo 



in 

■: I dlfg 




k 4 sin($ fc ) cos 2 (^fc) sin(?? 



g I e -2K kz za 



l^| 2 |r^)| 2 



(42) 



where 



and 



ki 2m e (e F -V ) ^ /In: \V„ 



r, 2 



(43) 



p (Kfc 2 ~iAe kq ) min(t,-t f ) / (ii kz ~iAe kq ) min(0, -t f ) 

r k3 (t) = Q(-t) — + e(t) - 



+ &{t- ti 



n kz - iAe kq \ K kz - iAe kq 

e -(«fc a +iAe fes )max(0,ty) _ e -(n kz +iAe kq )t \ 

' n kz + iAe kq J 

I 



(44) 



with 



Ae 



kg 



so 



£l - £k 



tq* 



zo 



8ir£ h 2 n q 2 cos 2 (i? 9 ) 



(45a) 
(45b) 



The term emerges from the inclusion of the transfer 
function © and needs to be set to zero when considering 
the population of the excited molecular level. 

Equation (|42|) is still highly complex, as it consists of 
an integral over 12 variables. In order to compute this 
expression numerically we divide the calculation into two 
steps. First, we calculate V~ on a discrete grid within the 

(k, <f)-space. Afterwards, we calculate the wave vector 
integral in Eq. (1421) while using sexa-linear interpolation 
to obtain the intergrid- values of V-jL. Because of this grid- 
based scheme we have to introduce a cut-off constant q c 



for the q- integration in Eq. (l4"2"j) . which we choose to be 
slightly larger than the classical cut-off q max - This is 
however not an issue, since the spectrum of the emitted 
electron falls off very rapidly beyond the classical cut- 
off. The high-dimensional integrals which occur in our 
formalism are then efficiently computable by means of 
Monte Carlo techniques. In particular, we employed the 
VEGAS algorithm as implemented in the GNU Scientific 
Library. 



V. RESULTS 

We now present numerical results. The parameter val- 
ues used in our calculations are listed in Table Q] and 
ILT1 In all of the calculations we fixed the turning point 
zq to the value given by Eq. @ for 50 meV, which 




FIG. 5: Final occupancies of the excited molecular level and 
of the free electron states (secondary electron emission co- 
efficient) in parallel (solid lines) and perpendicular (dashed 
lines) molecule orientation for different kinetic energies Skin 
of the incident molecule. The curves were calculated using 
the zeroth order formulas, Eq. ()23p and (|37a|) . respectively. 
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FIG. 6: Comparison of the secondary electron emission coef- 
ficient 7 e computed in zeroth order from Eq. (I37a|l (solid line) 
and in first order from Eq. (I37b|l (dashed line). The molecu- 
lar axis was aligned perpendicular to the surface. Clearly, the 
first order correction heals the divergence of 7 e for Skin — > 
and ensures 7e < 1 for all energies. 



amounts to approximately 4.42 as- Within the energy 
range < 1 eV this approximation is uncritical, since 
the turning point as given by Eq. Q varies very weakly. 
We aim at calculating, respectively, the occupancies of 
the excited molecular state n\ (t) and of the free electron 
states n(t) — J2q n q(t)- Within our model the initial 
magnetic quantum number fi = ±1 has no impact on the 
occupation numbers. For convenience we thus omit any 
m subscripts in the following. We restrict the molecule's 
orientation to the two fundamentally distinct situations 
(p = (axis parallel to the surface) and ip — ^ (axis 
perpendicular to the surface) . Furthermore, if not stated 
otherwise, we consider an aluminum surface. 

We start our analysis with the final occupancies, that 
is, rii(oo) and n(oo) where, according to Eq. (|34j) . the 
latter quantity is the secondary electron emission coeffi- 
cient 7 e . For now we employ only the lowest order equa- 
tions (|23|) and (|37a[) . which, concerning the emitted elec- 
tron, means that we neglect the matrix element renor- 
malization and thus the lifetime effects of the metastable 
molecule. 

The final occupancies of the excited molecular level 
and of the free electron states (secondary electron emis- 
sion coefficient 7 e ) are plotted in Fig. [5] for different ki- 
netic energies of the incident molecule. Obviously, the 
Penning process gets more efficient for lower kinetic ener- 



Parameter 


Value 


Reference 


S 


2.067 a B 


33 


eo 


-17.25 eV 


8 


El 


-9.57 eV 


1 



TABLE II: Model parameters for nitrogen. 



gies, which is evident, because lower kinetic energies cor- 
respond to smaller molecule velocities and thus to larger 
interaction times of the solid-molecule system. Note, that 
for Skin < I0 -6 eV the secondary electron emission coef- 
ficient gets larger than one although physical values for 
7e should obviously be less than or equal to one. This 
unphysical peculiarity is a consequence of the negligence 
of the lifetime effect in the zeroth order formula for j e , 
Eq. (I37a|) . and can be fixed by employing higher order 
terms of the full Eq. (|3l>|) (see below) . 

From Fig. [5] we see that for the two distinct orienta- 
tions the Penning process is almost equally efficient in de- 
exciting the molecule. The number of emitted electrons, 
however, is about a factor three smaller when the molec- 
ular axis is parallel to the surface as compared to when 
it is perpendicular. Numerical tests showed that for the 
parallel case the electron is primarily emitted with a very 
small perpendicular energy whereas in the perpendicular 
case the perpendicular energy of the emitted electron is 
distributed more equally. Thus, in the perpendicular ori- 
entation the electron has a higher probability to breach 
through the surface barrier originating from the image 
potential leading to a larger value for j e . 

To fix the unphysical behavior of the secondary elec- 
tron emission coefficient at low collision energies, the first 
order formula, Eq. (|37b|) . already suffices. It includes ad- 
ditional exponential factors which damp the integrand in 
the divergent region leading to "/ e — > f for ->■ 0. Fig- 
ure [6] explicitly demonstrates this behavior. The region 
in which the zeroth order result for 7 e exceeds unity cor- 
responds to very low molecule velocities Skin < 10~ 6 eV. 
These sub-thermal collision energies are rarely realized 
and are of lesser significance to our problem. In low- 
temperature plasmas, for instance, the systems we are 
primarily interested in, molecules have at least thermal 
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FIG. 7: Time evolution of the occupancies of the excited 
molecular level (upper panel) and of the free electron states 
(lower panel) for parallel (solid lines) and perpendicular 
(dashed lines) molecule orientation. The kinetic energy of the 
incident molecule was fixed to 50 meV . Time is measured in 
units of 1/(kv). A time difference of At = 1 thus corresponds 
to the motion of the molecule over a distance of Az — as/2. 
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FIG. 8: Energy spectrum n ef (oo) = n^(oo) of the Penning- 
emitted electron for parallel and perpendicular molecule ori- 
entation (reaction Q, this work). The kinetic energy of the 
incident molecule w as fixed to 50 meV. For comparison we 
also show the results fLorente et al.| £ obtained for reaction ((2}. 



collision energies. Only in beam experiments with ex- 
treme grazing incidence^ may the collision energies be 
low enough to require the lifetime effect to be explicitly 
included in a theoretical analysis. For electron energies 
above 0.1 meV the difference between the zeroth and 
first order formulas are vanishingly small indicating that 
in this energy range Eq. (|37aj) is sufficient. 

Next, we investigate the time evolution of the occupan- 
cies. We fix the kinetic energy of the molecule to 50 meV, 
which is about twice the thermal energy at room temper- 
ature. In addition, we employ the zeroth order formulas. 
This is justified because the higher order corrections are 
small in the considered energy region. The results are 
plotted in Fig. [7] Obviously, the occupancy of the ex- 
cited molecular level (upper panel) changes significantly 
only in the range \t\ < 5, which, taking the turning point 
into account, equals maximum distances of the molecule's 
center of mass from the surface of roughly 7 as ■ As ex- 
pected, within our model, the process is equally effective 
in the incoming and outgoing branch of the trajectory. 
The time evolution of the occupancy of the free states 
(lower panel) is however distinctively different. It shows 
a plateau around t = 0, that is, a stagnation of the prob- 
ability to escape from the solid. This is a consequence 
of the image potential which almost completely traps an 
electron emitted at low surface distances where its per- 
pendicular energy is too small to overcome the image 
barrier encoded in the surface transmission function. 

The energy distribution n e? of the emitted electron at 
t = oo is also of interest. This quantity is shown in Fig. [5] 
for the two principal molecule orientations and a collision 
energy Skm — 50 meV. The graphs for the two different 
orientations (denoted by "this work") start at the ori- 
gin and monotonously increase until a cut-off energy is 



reached. The latter resembles the classical cut-off energy 
£ max ^ gee pig P2]) ; implying that energy conservation is re- 
stored at the end of the collision, as it should be. The low 
energy part of the spectrum is cut off due to the surface 
transmission function which allows electrons to escape 
from the surface only when their perpendicular energy is 
large enough. The spectrum for the perpendicular case 
takes on larger values, as was explained in connection 
with the (^-integrated spectrum shown in Fig. [5] 

For comparison , we also plotted in Fig. [5] the spec- 
tra |Lorerrte^et7^yI obtained for an electron released due 
to charge-transfer and subsequent auto-detachment, re- 
action (|2]l. assuming zq s» 5a b- For low energies the 
spectra due to Penning de-excitation (JXJ) are practically 
zero and thus significantly smaller than the ones corre- 
sponding to reaction ([2]l, In this energy range the res- 
onant channel dominates. Close to the cut-off energy, 
however, the Auger spectra are only one order of mag- 
nitude smaller than the spectra for the charge-transfer 
process indicating that in this energy range the two pro- 
cesses are indeed eye-to-eye competitors. To study in our 
model the competition between the two channels quan- 
titatively requires however to dress-up the free Green's 
functions of the molecular electron by a hybridization 
self-energy and is thus beyond the scope of the present 
investigation. Integration of the energy spectra yields the 
secondary electron emission coefficients j e listed in Table 
IIIII For aluminum the Penning process is thus one to two 
orders of magnitude weaker than the resonant process. 

We now investigate the influence of the molecule's 
turning point on the de-excitation process. For alu- 
minum and tungsten the occupation number of the ex- 
cited molecular level calculated with the turning point 
given by Eq. (j4|) for s^in = 50 meV is shown in Fig. [9] 
together with the data obtained for the turning point 
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FIG. 9: Comparison of the final occupancy of the excited 
molecular level ni(oo) for two different values of the turn- 
ing point zo for aluminum (solid lines) and tungsten (dashed 
lines). The molecule's axis is perpendicular to the surface. 



fixed to zq — A as- For the latter the de-excitation prob- 
ability increases drastically at low energies leading to 
significantly reduced occupancies of the excited molec- 
ular state. At large energies the difference in the turning 
points shows almost no influence. Quantitatively, for en- 
ergies below 1 meV, the occupancy of the excited molec- 
ular level for zq — 4as is approximately one order of 
magnitude smaller as the one obtained for zq — 4.42 as ■ 
This has however no influence for the occupancy of the 
free states (not shown in Fig. [9]) because an electron emit- 
ted closer to the surface has a lower perpendicular energy 
and can thus no longer breach through the surface po- 
tential arising from the image potential as we discussed 
before. 

Throughout our investigation we assumed that the 
Coulomb interaction, which drives the Penning de- 
excitation of the molecule, is unscreened. In reality, how- 
ever, the Coulomb interaction in the vicinity of a sur- 
face is screened due to the charge carriers of the solid. 
The strength of Penning de-excitation should thus be af- 
fected by screening. To estimate this effect it suffices 
to consider the statically screened Coulomb potential 
^c( r ) = ^c( r )e _Ksr with k s the screening wave num- 
ber at the surface. Little is known about this quantity 





7e (parallel) 


7 e (perpendicular) 


This work 


1.93 • 1(T 3 


6.43 • 10 -3 


Lorcnte ct al. ' 


1.06 • 10" 1 


8.98 • 10~ 2 



TABLE III: Comparison of the secondary electron emission 
coefficient 7 e due to Penning de-excitation (this work) and 
resonant charg e-transfer with subsequent auto-detachment 
(jLorente et al.| £) for a A r 2( 3 S^") molecule hitting an aluminum 
surface with a kinetic energy of 50meV and parallel or per- 
pendicular orientation. 



except that it has to be smaller then the bulk screening 
wave number. Positron transmission and trapping exper- 
iments for various metallic films 36 indicated, for instance, 
that the screening wave number near the surface is most 
probably a factor 0.6 less than in the bulk. Taking this 
correction factor into account the screening wave num- 
ber for an aluminum surface, for instance, is 1.2285/A 
The screened Coulomb potential affects only the values 
of V£ ? , which enters quadratically into the lowest order 
formulas for the occupation numbers and reduces them 
by approximately 40%. Thus, as expected, screening re- 
duces the efficiency of the Coulomb-driven de-excitation 
channel, but it does not change its order of magnitude. 

So far, we co mpared our re sults only with the theo- 
retical results of lLorente et alj £ who studied the charge- 
transfer channel ([2]) for A/2( 3 E+) hitting an aluminum 
surface. To make conta ct to the experimental data ob- 
tained by IStracke et al.l 6 -. who also emphasized the di- 
rect charge-transfer reaction but nevertheless gave an 
estimate for 7 e due to Penning de-excitation ([1]), we 
also calculated j e for tungsten using the 50 meV turn- 
ing point of the surface potential (|4| and a statically 
screened Coulomb interaction with k s = 0.6 where 
= 1.7659/A is the screening wave number of bulk 
tungsten^ We find j e ss 2 • 10~ 3 which is rather close 
to the experimental estimate 7e Xp |ponnin g ~ 10~ 4 — 10~ 3 . 
Despite the crudeness of the effective model, which ne- 
glects, for instance, dangling bonds and surface states 
and the simplistic treatment of screening near a surface, 
our approach seems to capture even quantitatively the 
essential physics of Penning de-excitation. We attribute 
this to the rather large value of the molecule's turning 
point zq ~ 4.42as which partly immunizes Penning de- 
excitation against the details of the electronic structure 
in the immediate vicinity of the surface. 



VI. CONCLUSIONS 

We investigated the release of secondary electrons due 
to Auger de-excitation of metastable nitrogen molecules 
at metallic surfaces using an effective model for the two 
active electrons involved in the process and the Keldysh 
formalism to calculate the occupation numbers of the 
relevant single-electron states as originally proposed by 
iMakoshil ^ In contrast to him, however, we are not re- 
stricted to time-local Auger self-energies and thus to the 
wide-band approximation because we solve the Dyson 
equation for the retarded Green's function by exponen- 
tial resummation. We also employed the Auger matrix 
element obtained from a LCAO-type description of the 
molecule and an abrupt half-space- type description of the 
metal and not a phenomenological matrix element. The 
dependencies of the matrix element on the single-electron 
quantum numbers are retained in our calculation, as is 
the image interaction of the Auger electron with the sur- 
face and the distortion of the continuum of Auger states 
in the vicinity of the surface due to the molecule which 
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is accounted for by using for the Auger electron a con- 
tinuum of two-center Coulomb waves. 

The complexity of the final equations forced us to cal- 
culate and interpolate the Auger matrix element on a grid 
in a high-dimensional parameter space. This numerical 
approximation for the matrix element enabled us how- 
ever to calculate the time evolution and the final values 
of the occupancies of the excited molecular level and of 
the free electron states with standard Monte-Carlo inte- 
gration routines without further approximations. Since 
the Auger interaction is rather weak, we utilized only the 
lowest order formulas derived from the quan tum kine tic 
theory. The lifetime correction introduced bv lMakoshil ^ 
contained in higher order terms, was shown to be impor- 
tant only for very low kinetic energies of the molecule. 

We applied our approach to Auger de-excitation of 
S«) on aluminum or tungsten. For an aluminum 
surface we showed that for realistic turning points of the 
molecule's trajectory, for instance, the one obtained from 
the molecular-surface potential for thermal kinetic ener- 
gies, which is about 4.42as, Auger de-excitation is one to 
two orders of magnitude less efficient in releasing an elec- 
tron than the direct charge-transfer process investigated 
by Lorente and coworkers^ For a tungsten surface our 
model produced for a turning point of 4.42as a secondary 
electron emission coefficient due to Auger de-excitation 
which agrees with the experimental estimates of Strackc 
and coworkers. 6 

The effective model we used works with crude wave 
functions, but has the virtue to be parameterizable with 
a few easily obtainable energies. For the applications 
we have in mind, secondary electron emission in low- 
temperature gas discharges, where a great variety of dif- 
ferent kinds of molecules and different kinds of wall ma- 
terials occur, we consider this as a real advantage. With 
appropriate modifications the model can be applied to 
dielectric surfaces as well. We used the model only for 
the investigation of Auger de-excitation. In combina- 
tion with the Keldysh Green's function technique, how- 
ever, it can be also employed for the description of direct 
charge-transfer processes provided the energy shift and 
the broadening of the molecular levels due to the im- 
age interaction of the molecule with the surface and the 
self-energy corrections due to the Coulomb interaction 
between the excess electron and the excited electron of 
the molecule are included. 
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Appendix A: Wave functions 



The bound 



electron 



ive functions for the active molecular 
are calculated from the linear combi- 



nation of atomic orbitals 



^2p m (l)±^2 Pm (2) 



(Al) 



with 1 and 2 labeling the two distinct nitrogen atoms and 
N2-K n , g denoting normalization constants. The atomic 
nitrogen wave functions ^2p m are approximated by us- 
ing a hydrogen-like model with effective nucleus charge 
Z c ff = 4. The resulting wave functions show excellent 
agreement with Roothaan-Hartree-Fock calculations for 
the nitrogen atom.— The molecular wave functions (jAlj) 
are most conveniently expressed in cylindrical coordi- 
nates {g, tp, z) and possess the explicit form 



1 — 2m«2 



■ge 



(A2) 



x e 



-K^ + ( Z +f) 2 ± e - KV W(,-f) ; 



with k — 2 /as and normalization constants 



2±2k 5 / dg dz 

JO J-oo 



,-«vW(*+4) 2 P -<Ve 2 +(*-4) 2 



(A3) 



which need to be calculated numerically. The order of 
the wave functions within (IA1 1) . that is, the labeling of 
the nitrogen atoms does not matter in our calculation. 

The wave functions for the electrons inside the metal 
are calculated along the lines of Ref. [25| by solving the 
Schrodinger equation for an electron trapped by the step 
potential 



V(z) = 



-\V Q \ z<0 







z > 



(A4) 



The presence of the molecule is thus ignored, as far as 
the calculation of the metal wave functions is concerned. 
Using box normalization with box size L we obtain 

^r-(f) = —^ e i(Kx+k y y)( T ^ e -K kg z e ^ 



+ [e ik " + R k ,e^"] 0(-z)} , 

where the following wave vector dependent coefficients 
have been introduced 



ik z -\- K<k z 
ik z Kk z 

ik z Kk z 



2m e 



\Vo\~kl 



(A6a) 
(A6b) 
(A6c) 
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The energy of an electron in the state k is given by 

h 2 



e k 2m, 



(k 2 x +k 2 y +k 2 z )-\V Q \ 



(A7) 



We now turn to the free states involved in the Pen- 
ning process. As discussed in the main text we include 
here the effect the molecule has because the Auger elec- 
tron originates from the molecule. The continuum of free 
states resembles thus not the continuum of the solid but 
the single-electron continuum of the molecule which we 
approximate by a two-center Coulomb (TCC) wave^ In 
the past it has been successfully used to model electron- 
impact ionization and photo-ionization of H2 and 
molecules! 31 ' 32 The TCC wave is an approximate solu- 
tion of Schrodinger's equation for an unbound electron 
moving in the field of two fixed centers. Employing as 
as the unit of length it reads 



with 



(27T) 



NlC{q,r x ,Z l )NlC{q,ft l ,Z 2 ) , (A8) 



N% = e^i T[l+i— 



(A9a) 



C{q,r,Z)=M\-i-,l,-i[qr + q-r]\ , (A9b) 

where M is the confluent hypergeometric function of the 
first kind. The two vectors f\ and r*2 denote, respec- 
tively the position of the electron as seen from nucleus 
1 and 2 and Z\ and Z^ are effective charge numbers. As 
proposed in Ref. [3l| we account for the partial screen- 
ing of the nuclei by the passive electrons and choose 
Z\ = Z2 = Zc 



1 

2 ■ 



Appendix B: Keldysh formalism 

To fix our notation we give a brief description of the 
Keldysh formalism. For a more complete survey of the 
topic we refer the reader to Refs. [l4l - [l7l . 

Consider a fcrmionic system with a Hamiltonian 



H{t) = H + Hi{t). 



(Bl) 



where H\(t) represents a time-dependent perturbation 
of the non-interacting system Hq. Due to the time- 
dependence of the Hamiltonian H (t) we are faced with a 
non-equilibrium situation. 

One way to treat systems with a time-dependent 
Hamiltonian is the non-equilibrium Green's function 
technique introduced by KeldyshJ^ The key feature of 
the technique is a time contour C in the complex plane 
running from — 00 to 00 and then back again to —00. All 
quantities of the usual Green's function technique are 
then defined on this complex time path. 



Of particular importance is the contour-ordered 
Green's function G a p{t,t'), describing the propagation 
from a state (3 at time t' to a state a at time t. It is 
defined by 



Tr 



*a(*)*J(*' 



H 



(B2) 



where Tq specifies the chronological time-ordering oper- 
ator on the contour, and ^ represent the usual field 
operators, and (. . .)„ denotes the averaging with respect 
to an arbitrary state of the full dynamical system (|B1I) . 

Employing the interaction picture Eq. (|B2[) can be 
transformed to 



iQ a p{t,l/) = (Tch a (t)^)Sc\) 

\L J 'Ho 



(B3) 



with the contour scattering operator Sc defined by 

Sc =T c e-*Sc dt "^) . (B4) 



The tilde in Eq. (|B3|) and (|B4j) characterizes the cor- 
responding quantity in the interaction picture. Equa- 
tion (|B3[) is suitable for performing the usual perturba- 
tion expansion in terms of Hi(t), the only difference be- 
ing that all time integrals need to be taken over the time 
contour C instead of the real time axis. 

Inspecting the definition of the contour ordered 
Green's function (|B2[) and the possible location of the 
two time arguments t and t' on either the increasing (+) 
or the decreasing (— ) branch of the contour we can de- 
compose G a f}{t,t') into the four analytical pieces 



iG 



H 



a (5 (MO = 

iG-p(t,t') = (T a h a (t)*l(t' 



H 



H 



(B5a) 
(B5b) 
(B5c) 
(B5d) 



where T c and T a denote the chronological and anti- 
chronological time ordering operator on the real time 
axis, respectively. Equation (|B5j) can be also expressed 
in matrix notation, 



/■/-<++ r<+- 



G 



a0 



(B6) 



where the time arguments are omitted for convenience. 

The time evolution of the Green's function (|B6[) is gov- 
erned by the Dyson equation 



,(0) 



y aS ^Sj y iP ' 



(B7) 



where the summation over internal indices and integra- 
tion over internal times is implicitly assumed. Equa- 
tion (|B7[) involves the unperturbed Green's function (in- 
dicated by a (0) superscript) and the self-energy E<$ 7 . 



14 



The latter is defined on the contour C as well and thus 
possesses a matrix representation similar to Eq. (1B6[) . 
Hence, the Dyson equation (|B7j) is a matrix equation. 

The structure of this equation can be simplified when it 
is noted that the set (|B5I) is linearly dependent. Applying 
the unitary transformation 



U 



J_/l -1 
V^l 1 1 



(B8) 



that determines the different parts of the Green's func- 
tion, 



G 



A/R 
a/3 
K 



G 



a/3 



A/R(0) 

a/3 

K(0) 



G 



A/R(0)^A/R qA/R 



7/3 



a/3 

G 



n K(0)^A n A 
<^„x ^S-y 7/3 



R(0) 



x6 

^57 °7/3 



^R n K 

^7 ^7/3 



(Blla) 
(Bllb) 



to the Dyson equation (|B7[) the Green's function and the 
self-energy turn into 



Gaf3 = UQ a p J7 t : 
^a/3 — U Y, a p W 






r~iA 
u a/3 


r<R 


^a/3 


'yK 

^a/3 


yR 

^a/3 


Pa/3 






(B9a) 



(B9b) 



The superscripts A, R and K denote, respectively, the ad- 
vanced, retarded and Keldysh part of the corresponding 
quantity. For the retarded and advanced Green's func- 
tions holds 



Equation (|Bllb|) can be solved iteratively to give the 
important relation 



n K _ 

U al3 - 



a7 



n R v fl 



r K(0) 



G t/3 + <^/3 



r R yK n A 
U a(5 ^57 U 7/3 



(B12) 



The Keldysh part of the Green's function can thus be 
computed from Eq. (|B12|) once the advanced and re- 
tarded parts are known. The diagonal component G^ a 
can then be used to calculate the occupation of the state 
a at arbitrary times 



(BIO) 



Carrying out the matrix multiplication in the trans- 
formed Dyson equation one obtains the set of equations 



*a(*) = ^[l-iG* (*,*)] • 



(B13) 



* Contact: marbach@physik.uni-greifswald.de 
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